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Abstract. In this paper we study the phenomenon of nonUnear supratrans- 
mission in a semi-infinite discrete chain of coupled oscillators described by 
modified sine-Gordon equations with constant external and internal damping, 
and subject to harmonic external driving at the end. We develop a consistent 
and conditionally stable finite-difference scheme in order to analyze the effect 
of damping in the amount of energy injected in the chain of oscillators; nu- 
merical bifurcation analyses to determine the dependence of the amplitude at 
which supratransmission first occurs with respect to the frequency of the driv- 
ing oscillator are carried out in order to show the consequences of damping on 
harmonic phonon quenching and the delay of appearance of critical amplitude. 



1. Introduction 

The study of nonlinear continuous media described by modified Klein-Gordon 
equations subject to initial conditions is a topic of interest in several branches of 
the physical sciences [U [2j |3j H] . Here the analytical study on features of solutions 
of Klein-Gordon-like equations as well as the development of new and more pow- 
erful computational techniques to approximate them have been the most transited 
highways in the mathematical research of the area. 

The behavior of nonlinear continuous media described by modified Klein-Gordon 
equations subject to boundary conditions is also an interesting topic of study in 
mathematical physics. This type of mathematical models have proved to be useful 
when applied, for instance, to the description of the project of data transmission 
in optical fibers in nonlinear Kerr regimes [5j [6] or in the study of the property of 
self-induced transparency of a two-level system submitted to a high-energy incident 
(resonant) laser pulse [3 [8]. More concretely, the behavior of a continuous medium 
submitted to a continuous wave radiation constitutes a fundamental problem that 
has not been studied in-depth, yet it has shown to have potential applications as a 
model in the study of Josephson transmission lines [HI UHl E] • 

Numerical studies on discrete versions of these models [12] followed by experi- 
mental results 13 have shown that there exists a bifurcation of wave transmission 
within a forbidden band gap in certain semi-infinite undamped Klein-Gordon-like 
chains of coupled oscillators which are periodically forced at the end. This bifurca- 
tion is a consequence of the generation of nonlinear modes by the periodic forcing 
at the end of the chain, and allows energy to flow into the medium via a nonlinear 
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process called nonlinear supratransmission, which has been proved numerically not 
to depend on integrability as long as the model possesses a natural forbidden band 
gap. 

The process of nonlinear supratransmission has been studied numerically by 
means of computational algorithms already built in standard mathematical pack- 
ages. Most particularly, the numerical results obtained in '12', for instance, rely on 
the use of a Runge-Kutta method of high order, which has the advantage of possess- 
ing high accuracy and efficiency in the computations, but lacks the consistency in 
the numerical estimation of the energy flowing into the medium which is desired in 
the study of supratransmission. With this shortcoming in mind, we present in this 
paper an alternate numerical formulation to study the process of supratransmission 
in a nonlinear system of differential equations, that has the advantage of being con- 
sistent in the approximation of the solutions to the problem and in the estimation 
of the continuous energy. More concretely, in the present work we study the pro- 
cess of nonlinear supratransmission in a semi-infinite nonlinear discrete system of 
coupled oscillators governed by modified sine-Gordon and Klein-Gordon equations 
with constant internal and external damping. We develop some numerical results 
for dissipative nonlinear modified Klein-Gordon-like equations, and validate our 
conclusions against [12, . Our study — rigorous in numerical nature — will soon be 
succeeded by future applications in forthcoming papers. 

In Section II of this paper we introduce the mathematical problem under study 
and derive the expression of the instantaneous rate of change of the energy trans- 
mitted to the system through the boundary. Section III is devoted to introducing 
the finite-difference scheme; here we present the discrete energy analysis associated 
with the problem under study and establish in detail the numerical properties of 
our method. Numerical results are presented in the next section, followed by a brief 
discussion. 



2.1. Mathematical model. In this article we study the effects of the nonnegative 
constant parameters (3 and 7, and the real constant on the behavior of solutions 
to the mixed- value problem with mass m, 



which describes a system of coupled oscillators with coupling coefficient c » 1, and 
where /3 and 7 clearly play the roles of internal and external damping coefficients, 
respectively. The number A^it„ is used to represent the spatial second-difference 
Un^i—2u„+u„-i for every n € Z+, and the boundary-driving function ip is assumed 
continuously differentiable on (0,-|-cx)). In our study, we will let V{u) = 1 — cosu 
for a chain of coupled sine-Gordon equations, and V{un) = — j\U^ + for 
zero mass in the case of a Klein- Gordon chain. 

It is worth recalling that the differential equation under study has multiple ap- 
plications in the continuous limit. For instance, a similar equation appears in the 
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study of Josephson junctions between superconductors when dissipative effects are 
taken into account [1]. A modification of this equation is used in the study of 
fluxons in Josephson transmission hnes [2" , and a version in the form of a modified 
Klein-Gordon equation appears in the statistical mechanics of nonlinear coherent 
structures such as solitary waves (see [3^ pp- 298-309) when no internal damping 
is present. Besides, our equation clearly describes the motion of a damped string 
in a non-Hookean medium. 

For purposes of this paper, we will consider a system of coupled oscillators ini- 
tially at rest at the origin of the system of reference, subject to an external harmonic 
forcing described by "tpit) — Asinflt, and pure-imaginary or real mass satisfying 
|m| <^ 1. Analysis of the undamped linearized system of differential equations in 
^ shows that the linear dispersion relation reads a;^(fc) = + 1 + 2c^(l — cosfc) 
in any case. The driving frequency will take on values in the forbidden band gap 
region < vn? + l, in which case the linear analysis yields the exact solutions 



Unit) = 74sin(fit)e"^", where 



A = arccosh 1 



2c2 

Meanwhile, the massless undamped nonlinear case possesses an exact solution 
in the continuous limit which happens to work well for high values of the cou- 
pling coefficient. It has been shown numerically that, for each frequency f2 in the 
forbidden band gap, the massless medium starts to transmit energy by means of 
nonlinear mode generation for amplitudes greater than the threshold value As — 
4 arctan (Ac/fi). 

2.2. Energy analysis. In this section we derive the expression of the instantaneous 
rate of change of total energy in system ([1]) . Here, by a square- summahle sequence 
we understand a real sequence (an)^Q for which ^ is convergent. 

Lemma 1 (Discrete Green's First Identity). // (an)5^o square- summahle se- 
quence then 

oo oo 

^(a„+i - 2a„ + a„„i)a„ = ao(ao - ai) - ^(a„ - a„_i)^. 

n— 1 n—1 

Proof. Holder's inequality implies that both series in the formula of the lemma 
converge. Moreover, the sequence defined by t„ = a„a„_i — (a„_i)^ for every 
positive integer n converges to zero and the result follows now from the identities 



^^(On+l — 2a„ + an-l)a„ — y~^(^Tt+l — tn) — y~^(aK — On-l)^ 

= 1 n=l 

oo 

tl - ^(fln - an-l f 



71—1 n—1 n—1 

^2 



□ 



Proposition 2. Let (un(i))^o solutions to ^ such that {un{t))^^Q is square- 
summable at any fixed time t. The instantaneous rate of change of the total energy 
in the system is given by 



— = c^(uo - ui)uo ~ 7y^("«) 

n—1 



oo 

E 

,n=l 



(Un - Un-l)^ + (Ul - U0)U0 
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Proof. The energy density of the undamped system of coupled equations with a 
potential energy V{un) in the n-th oscillator is given by iJ„ = + c^(u„+i — 
w„)^ + rn^u^] + V{un)- After including the potential energy from the coupling 
between the first two oscillators, the total energy of the system at any time becomes 

oo 2 

and the fact that u„ tends to as ri increases implies that the sequence J„ = 
—c^Un(un — Un-i) converges to zero pointwisely at any fixed time. Simplifying and 
rearranging terms of the derivative of the Hamiltonian iJ„ yields the expression 

dH 

= {Jn - Jn + l) + - 2u„ + U„_i)u„ - 7(u„)^. 

The result follows now from this identity after computing the derivative of the 
total energy of the system, using the formula for telescoping series and applying 
our discrete version of Green's First Identity. □ 

Observe from the proposition that the expression of the rate of change of energy 
associated with damped system ([1]) is in agreement with the undamped formula de- 
rived in [H]. Moreover, it is clear that the external damping coefficient contributes 
decreasing the total amount of energy in the chain system for /3 equal to zero. 



3. Numerical analysis 

3.1. Finite-difFerence scheme. From a practical point of view, we will restrict 
our study to a system consisting of a finite number N of coupled oscillators with 
constant internal and external damping coefficients, described by the system of 
ordinary differential equations 



d<2 



■/3 



dt 



AX 



du„ 
dt 



+ m^Uri + V'{Un) = 



for 1 < n < N , where a includes both the effect of external damping and a simula- 
tion of an absorbing boundary slowly increasing in magnitude on the last N — No 
oscillators. More concretely, we let UN+i{t) be equal to zero at all time t, and let 
a be the sum of external damping and the function 



7'W - 



1 + tanh 



2n- No + N 
2^ 



0, 



No <n< N, 
otherwise. 



In practice, we will let k = 0.5, a = 3, Nq = 150 and — 200. 

We proceed now to discretize problem ^ using a finite system of N differential 
equations and a regular partition = to < < • ■ ■ < = 7" of the time interval 
[0,T] with time step equal to At. For each k = 0, 1, . . . ,M, let us represent the 
approximate solution to our problem on the n-th oscillator at time tk by u^. If 



we convey that Stu„ 



,k+l 



,fc-l 



that S. 



'2„,k 



,k+l 



,k-l 



and that 
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(5^M„ — — 2m„ + u„_i, our problem takes then the discrete form 



(2) 



A. 

2 At 



m 



,k+l 



subject to 



,fc+i 



fc-i 



u\ = (l>{n 

k 



0, 



At(p(r 



Af+l 



= 0, 



1 < n < iV, 
1 < n < iV, 
l<k<M, 
l<k<M. 



Observe that the proposed numerical method is nonlinear and requires an ap- 
plication of Newton's method for systems of equations in order to be implemented. 
Notice also that if V'{un) at the fc-th time step is approximated by V'{u'^) then 
the finite-difference scheme becomes linear and an application of Grout's reduction 
technique for tridiagonal systems suffices to approximate solutions of ([T|). In such 
case, it is readily seen that the vector equation Au'^'^^ = Bu'^ + Cu^~^ — V {u'^)+\''' 
must be satisfied for every k > 2, for the three iV-dimensional tridiagonal matrices 
and the A^-dimensional vector 



C = 



( b a ■■■ 

a b 



V ••■ 

/ e a ■ ■ ■ 

a e 

V 







b J 







B = 



( d c2 • • • \ 
C2 d ■■■ * 



V 



dj 



( (?u%- aStU^ \ 





respectively, and constants 



d^ 



_J_ 
2 At' 
2 

jAty 



2c' 



and 



b = 



a + 2p 

2At 

a + 2p 

2At 



m 

9 

m 

'Y 



1 



(Ai)2' 

I 



Here u*^ = (mJ^, . . . , uj^)* for every k e {0, 1, . . . , Af}, and V{u^) is the n-th dimen- 
sional vector whose i-th component is equal to V{u\). This latter formulation of 
our problem will be used for validation purposes only, since we will prefer the non- 
linear formulation due to the quadratic order of convergence of Newton's method 
and other reasons that will be presented in the next section. 
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Figure 1. Approximate solution ueo{t) of a sine-Gordon system 
in the first row and a Klein-Gordon system in the second, for a 
driving frequency of 0.9 and two different amplitude values A: (a) 
before and (b) after the bifurcation threshold. 

3.2. Discrete energy. The total energy in the system at the k-th time step will 
be approximated numerically via 

11=1 ^ ^ n=l 

2 ^ 2 ^ 2 

n— 1 71—1 

Proposition 3. T/ie following identity holds for every sequence (ujij) satisfying 
scheme ^ and every positive index k: 

At ~ "i^2Ai ^^l2At 

n=l ^ 



E 



2Ai y V 2At j 2Ai 

.n— 1 
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Figure 2. Graphs of total energy transmitted into the massless 
sine-Gordon and Klein-Gordon systems vs. driving amplitude for 
a driving frequency of 0.9, with P — and 7 = (sohd), 0.01 
(dashed), 0.02 (dash-dotted), 0.03 (dotted). 



Proof. The proof of this result is merely an algebraic task. We need only say that 
an application of the discrete Green's First Identity with a„ = {u'^'^^ — 2u^~^)/2Ai 
is indispensable in order to reach the correct expression of the term with coefficient 
p. □ 

We conclude that the discrete energy associated with scheme (H]) is a consistent 
approximation of order 0{At) of the total energy of system ([T]), whereas the discrete 
rate of change of energy is consistent order 0{At)^ with the instantaneous rate of 
change. 

3.3. Stability analysis. The following result summarizes the main numerical prop- 
erties of our method. 

Proposition 4. Scheme ^ is consistent order 0{At)^ with the linear contribution 
of ([TJ for a potential equal to zero. Moreover, a necessary condition for the scheme 
to be stable order n is that 

c^-^) (At)^<l+(f + /3)A.. 

Proof Following the notation in [14j, let U^^^ = 1*^;+^ and U2n^ = u'^ for each n = 
0,1,..., M and k = 0,1, . . . ,N - 1. For every n = 0, 1, . . . , M and fc = 1, 2, . . . , 
let be the column vector whose components are and [/^n- Our problem can 
be written then in matrix form as 

where 



1/" \ 1 ' 



At ^At^2 2(^0 

q = 1 + a B — Si+m^- — and 

^ 2 2 2 

At ^At^2 2(^0^ 

h = 1-a h/? — Si+m^- — —. 

2 2 2 
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Figure 3. Total energy transmitted into the sine-Gordon system 
vs. driving amplitude for a driving frequency of 0.9, with /3 = 0, 
7 = 0.01, and pure-imaginary and real masses of magnitude 
(solid), 0.05 (dashed), 0.075 (dash-dotted) and 0.1 (dotted). 



Applying Fourier transform to the vector equation we obtain 
U^^^=( |§y(l-2cWsin^|) -fi 

The matrix A(^) multiplying in this equation is the amplification matrix of 
the problem. It is easy to check that the eigenvalues of A for £, — tt are given by 

1 - 2c2(At)2 ± y'(l-2c2(Ai)2)2_^(7r)5(7r) 



Suppose for a moment that 1 — 2c^(Af)^ < — ^(Tr). If the radical in the expression 
for the eigenvalues of A{t:) is a pure real number then | A_ | > 1. So for every n g N, 
1 1 ^"11 > I'^-l" grows faster than Ki +nK2 for any constants Ki and K2- A similar 
situation happens when the radical is a pure imaginary number, except that in 
this case | • | represents the usual Euclidean norm in the field of complex numbers. 
Therefore in order for our numeric method to be stable order n it is necessary that 
1 — 2c^(At)^ > — ^(tt), which is what we wished to establish. □ 



4. Numerical results 

In this section we study the effects of the internal and external damping coef- 
ficients on the behavior of solutions to mixed- value problem ([T]). Particularly, we 
wish to establish the effect of weak damping on the minimum amplitude value nec- 
essary for the phenomenon of supratransmission to take place at a fixed driving 
frequency. Throughout we validate our code against [12, and against an imple- 
mentation of our finite-difference scheme using the Runge-Kutta method of order 
four. 

4.1. External Damping. For the remainder of this work we consider a semi- 
infinite coupled chain of oscillators initially at rest in their equilibrium positions, 
subject to harmonic forcing at the end described by ^{t) — Asin{V,t) at any time 
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Figure 4. Graphs of total energy transmitted into the massless 
sine-Gordon and Klein-Gordon systems vs. driving amplitude for 
a driving frequency of 0.9, with 7 = and (3 — (solid), 0.1 
(dashed), 0.2 (dash-dotted), 0.3 (dotted). 



t. The functions (j) and (f are both identically equal to zero and, in order to avoid 
the shock wave produced by the vanishing initial velocity in the boundary, we let 
the driving amplitude linearly increase from to A in a finite period of time before 
we start to compute the total energy. In the present section we will let (3 be equal 
to zero and consider a discrete system of 200 coupled oscillators described by (P) 
over a typical time period of 200, with a time step of 0.05 and a coupling coefficient 
equal to 4. 

Let us first examine the massless case when no damping is present at all and 
f2 — 0.9. In order to verify that our mixed-value problem produces a bifurcation 
it is necessary to study the qualitative behavior of solutions near the predicted 
threshold value As which, in this case, reads approximately 1.80. The first row 
of Figure [T] shows the function ugQ{t) in the solution of a sine-Gordon system for 
two different values of the amplitude of the driving end, while the second row 
presents the corresponding graphs in the solution of a Klein-Gordon chain. The 
graphs evidence the presence of a bifurcation in the behavior of solutions around 
the amplitude value 1.79 for both chains. 

Naturally, the next step in our investigation will be to determine the behavior of 
the total energy flow injected by the periodic forcing at the end of the undamped 
discrete chain of oscillators as a function of the amplitude. The solid lines of Figure 
[5] represent the graphs of total energy transmitted into the system vs. amplitude 
for a forcing frequency of 0.9 and external damping coefficient equal to zero, for a 
sine-Gordon system in the first column and a Klein-Gordon system in the second. 
It is worthwhile noticing the abrupt increase in the total energy administered to the 
system for some amplitude value between 1.77 and 1.79. This feature of the graphs 
evidences the existence of a bifurcation value around the predicted amplitude A^, 
after which the phenomenon of nonlinear supratransmission takes place. 

Figure [2] also presents graphs of total energy vs. forcing amplitude for weak 
constant damping coefficients 7 = 0.01, 0.02 and 0.03 in a sine-Gordon system of 
oscillators. The graphs show a tendency of the bifurcation value to increase linearly 
as the external damping coefficient is increased. Another interesting feature of this 
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Figure 5. Total energy transmitted into the sine-Gordon system 
vs. driving amplitude for a driving frequency of 0.9, with /3 = 0.2, 
7 = 0, and pure- imaginary and real masses of magnitude (solid), 
0.05 (dashed), 0.075 (dash-dotted) and 0.1 (dotted). 



figure is the decrease of total energy for increasing values of 7, at least for fixed 
amplitudes greater than the undamped bifurcation threshold. Needless to point 
out that similar conclusions are obtained for Klein-Gordon chains of oscillators. 

The qualitative effect of m in a sine-Gordon system is also of interest in the 
analysis of solutions of this chain and is numerically carried out in Figure [3] for 
a chain with external damping equal to 0.01 and pure-imaginary and real masses, 
using graphs of total energy administered into the system through the boundary vs. 
driving amplitude. A horizontal shift in the occurrence of the bifurcation value is 
readily noticed in these graphs. Indeed, the displacement of the bifurcation ampli- 
tude for the system ([T|) of mass m with respect to the bifurcation amplitude of the 
massless system is a monotone increasing function of m? . Analogous computational 
results (not included here) were obtained for a similar Klein-Gordon chain. 



4.2. Internal damping. Consider again a system of 200 oscillators ruled by mixed- 
value problem ([T]) over a time period of 200, with time step 0.05, coupling coefficient 
equal to 4 and constant external damping equal to zero. In this context. Figure |4] 
shows graphs of total energy vs. forcing amplitude for internal damping coefficients 
/? — 0, 0.1, 0.2 and 0.3, for sine- Gordon and Klein- Gordon systems. As in the case 
of external damping, we observe that the threshold value at which supratransmis- 
sion starts tends to increase as the value of /3 is increased. Opposite to the case of 
external damping, though, the minimum value for which supratransmission starts 
varies with /3 in a nonlinear way. 

Next we verify our conclusions on the effect of the mass m on the qualitative 
behavior of solutions of the sine-Gordon system. Figure [5] shows graphs of total 
energy vs. driving amplitude for a system with no external damping, /3 = 0.2 and 
driving frequency of 0.9, for different pure-imaginary and real masses. As observed 
before, the graphs evidence a shift on the bifurcation amplitude with respect to the 
massless bifurcation value which is an increasing function of m? . 
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Figure 6. Graphs of total energy vs. driving frequency and am- 
plitude for massless and undamped sine-Gordon and Klein-Gordon 
chains of coupled oscillators. 



4.3. Bifurcation analysis. Both the chain of coupled sine-Gordon equations and 
the chain of Klein- Gordon equations have numerically proved to undergo nonlin- 
ear supratransmission for a driving frequency equal to 0.9 and different values of 
the external and internal damping coefficients, thus establishing that the results 
obtained in this work do not depend on integrability. Naturally we are interested 
in determining that the process of nonlinear supratransmission happens for any 
frequency value in the forbidden band gap. With that purpose in mind, we ob- 
tained graphs of total energy vs. driving frequency and various amplitude values 
for undamped discrete systems of 200 coupled oscillators with coupling coefficient 
4, time period of 200 and time step 0.05. The 3-dimensional results for both chains 
of oscillators are shown in Figure [S] together with a graph of the continuous-limit 
threshold amplitude As vs. driving frequency on the amplitude-frequency plane for 
comparison purposes. 

Several observations may be immediately drawn from Figure [51 First of all, 
the predicted bifurcation values As display an excellent concordance to the cor- 
responding values obtained using the finite difference-schemes associated with the 
sine-Gordon and the Klein-Gordon chains. Second, the process of supratransmis- 
sion ceases to appear in the Klein-Gordon case for driving frequencies below 0.7. 
Third, for driving frequencies close to the band gap limit, the bifurcation thresh- 
old is not clearly determined from the energy vs. driving amplitude graph of the 
sine-Gordon chain, as prescribed by [12]. For those frequencies, it is indispensable 
to increase the time period of approximation at least up to 500. Fourth, strong nu- 
merical proof of the existence of the occurrence of the supratransmission process is 
at hand and we have now reasons to believe that there exists a bifurcation function 
A(ri;/3,7,m^) for driving amplitude associated with ([T]). 

We proceed then to obtain graphs of amplitude values for which nonlinear supra- 
trasmission starts vs. driving frequency for a massless sine-Gordon chain of coupled 
oscillators with internal damping coefficient equal to zero and different values of 
7. The numerical results are summarized in Figure [7] together with the plot of the 
prescribed continuous-limit bifurcation amplitudes As. It is worth noticing that the 
bifurcation threshold increases with 7 for fixed frequencies above 0.35, as previously 
evidenced for a frequency equal to 0.9. We must notice also that the discrepancies 
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Figure 7. Diagram of smallest driving amplitude value at which 
supratransmission begins vs. driving frequency for a massless sys- 
tem with /3 = and values of 7 equal to (solid), 0.1 (dashed), 
0.2 (dash-dotted) and 0.3 (dotted). The theoretical threshold As 
in the undamped is shown as a sequence of x -marks. 

that appear for frequencies below 0.35 when 7 = 0, also appear for greater values 
of 7, each time bounded in smaller intervals. In fact, we have checked that the 
discrepancy region — an effect of the phenomenon of harmonic phonon quenching 
— tends to vanish for higher values of external damping (results not included). In 
this state of matters, we wish to point out that better numerical approximations 
to the bifurcation threshold are obtained for larger systems of oscillators at the 
expense of superior needs in terms of computational time. Likewise, we have con- 
firmed that larger values of the coupling coefficient lead to better approximations 
to the bifurcation threshold, as prescribed by |12) . 

Figure [5] shows bifurcation diagrams for a massless sine- Gordon system of oscilla- 
tors with no external damping and different constant internal damping coefficients. 
As expected, the bifurcation threshold tends to increase with /? for a fixed driving 
frequency. The effects of harmonic phonon quenching are present again in all the 
bifurcation diagrams, but contrary to the case of external damping, in the case 
of internal damping the range over which discrepancies occur slightly widens as /3 
increases. Also, it is worth observing that the length of the forbidden band gap 
increases with the parameter /? apparently in a linear fashion. Moreover, the graphs 
of bifurcation diagrams for nonzero values of (3 are approximately obtained by shift- 
ing horizontally the corresponding graph of the undamped diagram a number of 
/3 units to the right. More concretely, A(ri; 0,0,0) is a good approximation for 
A(r!-/3;/3,0,0). 

Finally, we compute bifurcation diagrams for an undamped sine-Gordon system 
with different real and pure-imaginary masses in order to establish the effect of m 
on the occurrence of the bifurcation threshold. The numerical results are summa- 
rized in Figure [9] for some real and pure- imaginary masses, and driving frequencies 
starting at 0.5. Our results lead us to conclude that the bifurcation diagram of 
the sine-Gordon system of oscillators with mass m is approximately equal to the 
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Figure 8. Diagram of smallest driving amplitude value at which 
supratransmission begins vs. driving frequency for a massless sys- 
tem with 7 = and values of /? equal to (solid), 0.1 (dashed), 
0.4 (dash-dotted) and 0.7 (dotted). The theoretical threshold As 
in the undamped case is shown as a sequence of x -marks. 



graph of the corresponding massless graph shifted ^/^ + rn? — 1 horizontal units, 
for |m| ^ 1. In order to test our claim numerically, we obtained graphs of abso- 
lute differences between the massless undamped bifurcation diagram, and shifted 
undamped bifurcation diagrams for several mass values (results not included). We 
observed that the differences tend to attenuate for high frequency values and that 
smaller differences are obtained for smaller values of m in magnitude. 



In this paper we have developed a numerical method to approximate solutions 
of a semi-infinite nonlinear chain of coupled oscillators ruled by modified sine- 
Gordon equations harmonically driven at its end. The proposed finite-difference 
scheme is consistent order C(At)^ and we provided a necessary condition in order 
for the method to be stable order n. The process of nonlinear supratransmission 
for a coupled system of oscillators described by sine-Gordon equations was studied 
numerically under the scope of this numerical technique, and the dependence of 
supratransmission on damping was analyzed. 

Several conclusions can be drawn from our computational experiments on the 
sine- Gordon system of coupled oscillators ([T]). First of all, we have shown that the 
phenomenon of harmonic phonon quenching still appears in the presence of external 
and internal damping and that the discrepancy region due to phonon quenching is 
shortened as the external damping coefficient is increased, while it slightly widens 
as the internal damping coefficient increases. Second, the threshold value at which 
supratransmission first occurs for fixed frequencies outside the discrepancy region is 
seen to increase for both external and internal damping as the damping coefficient 
increases; both conclusions are clearly consequences of the dispersive and dissipative 
natures of the parameters and 7, respectively. 
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Figure 9. Diagram of smallest driving amplitude value at which 
supratransmission begins vs. driving frequency for an undamped 
sinc-Gordon system with -\/m| + 1 = 1 + for £ = —4, —2, — 1, 
0, 1,2,4. 



Finally, the bifurcation diagram for a value of the parameter fH is approximately 
equal to the corresponding diagram for the undamped system shifted /3 horizontal 
units to the right. Likewise, a horizontal shift of ^/l + w? — ! units in the bifurcation 
diagram of a sine-Gordon system of mass m with respect to the corresponding 
massless system is observed for small masses and frequencies outside the discrepancy 
region. 
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